2203.10007v1 [astro-ph.SR] 18 Mar 2022 


. 
. 


arXiv 


OB Associations 


Nicholas J. Wright 
Astrophysics Group, Keele University, Keele ST5 5BG, UK 


Simon Goodwin 
Department of Physics and Astronomy, University of Sheffield, Sheffield, S3 7RH, UK 


Robin D. Jeffries 
Astrophysics Group, Keele University, Keele ST5 5BG, UK 


Marina Kounkel 


Department of Physics and Astronomy, Vanderbilt University, VU Station 1807, Nashville, TN 37235, USA 


Eleonora Zari 
Max-Planck-Institut für Astronomie, Königstuhl 17 D-69117 Heidelberg, Germany 


OB associations are low-density groups of young stars that are dispersing from their birth 
environment into the Galactic field. They are important for understanding the star formation pro- 
cess, early stellar evolution, the properties and distribution of young stars and the processes by 
which young stellar groups disperse. Recent observations, particularly from Gaia, have shown 
that associations are highly complex, with a high degree of spatial, kinematic and temporal 
substructure. The kinematics of associations have shown them to be globally unbound and ex- 
panding, with the majority of recent studies revealing evidence for clear expansion patterns in 
the association subgroups, suggesting the subgroups were more compact in the past. This expan- 
sion is often non-isotropic, arguing against a simple explosive expansion, as predicted by some 
models of residual gas expulsion. The star formation histories of associations are often complex, 
exhibit moderate age spreads and temporal substructure, but so far have failed to reveal simple 
patterns of star formation propagation (e.g., triggering). These results have challenged the his- 
torical paradigm of the origin of associations as the expanded remnants of dense star clusters 
and suggests instead that they originate as highly substructured systems without a linear star 
formation history, but with multiple clumps of stars that have since expanded and begun to over- 
lap, producing the complex systems we observe today. This has wide-ranging consequences for 
the early formation environments of most stars and planetary systems, including our own Solar 


System. 


1. Introduction 


OB associations are low-density groups of young stars, 
typically containing many prominent OB stars as well 
as numerous low-mass stars 
[1999). They are technically divided into OB associations, 
which contain bright OB stars, and T associations, contain- 
ing prominent T-Tauri stars, though other than their total 
mass there are no other differences and the term associa- 
tion is often used to refer to both types. 

The low space densities (< 0.1 Mo pc?) of associ- 
ations make them dynamically unstable to Galactic tidal 
forces and therefore over time they should disperse. The 
fact that they exhibit some spatial and kinematic concentra- 
tion (despite most likely being unbound) and contain short- 
lived OB stars implies that they must be young, and are 
therefore valuable tracers of the star formation process, al- 


lowing us to study the propagation of star formation and the 

role of feedback in triggering or halting star formation. 
Historically, OB associations were vital for identifying 

groups of OB stars and calibrating their luminosity scale 


1953 1978). This work led to 


the first census and catalogue of classical OB associations 
by [Ruprecht (1966), which was compiled and standardised 
from earlier works. This catalogue was most-recently up- 
dated by [Wright] (2020), though many systems still remain 
poorly defined and studied. More recently, associations 
have provided large samples of unobscured young stars that 
are useful for studies of the initial mass function, the fre- 
quency and properties of multiple systems, protoplanetary 


disks and planetary systems (Massey et al.|1995 
hoven er al [2007] Kalas er al. 2013). 


The origin of associations is still debated; according to 


some they are the short-lived expanded remnants of dense 
star clusters (Kroupa|2011), while others suggest they are 
the result of multiple, sub-structured star formation events 
occurring over large spatial scales or at low-density 


and Scalo[I978] [Krujjssen]2012). 


1.1. Overview of association properties 


Associations have total stellar masses between a few 
hundred and a few tens of thousand of solar masses. Their 
dimensions range from a few tens to a hundred parsecs 
(Section 3.1), giving them typical densities of 0.001 — 0.1 
Mo pe’. They are highly asymmetric and substructured 
(Section [.1], and often contain open clusters or star form- 
ing regions within their boundaries. Associations have 1D 
velocity dispersions of a few km s^! and are generally ex- 
pected to be gravitationally unbound (Section 3.2) and ex- 
panding (Section[3.3p. 

The ages of associations range from a few to a few tens 
of Myr, with both considerable age spreads (comparable to 
the age of the system) and resolvable age substructure (Sec- 
tions[3.4]andJ3.5p. The lower age limit historically separated 
associations from obscured and embedded systems (embed- 
ded clusters, [Lada and Lada|2003), while the upper limit 
comes from the difficulty identifying low-density groups of 
stars older than this (Section[2.1) — though see[Kounkel and] 
for an attempt to extend this to older systems. 


Associations are observed throughout the disk of our 
galaxy and their distribution has been used to trace the spi- 
ral structure of the Milky Way (Section 3.6), map previous 
generations of young stars and study the propagation of star 
formation (Section[4.2). 


1.2. On the definition of OB associations 


OB associations were originally defined as systems of 
stars that are very young compared to the age of the galaxy, 
have a common origin, and a stellar density lower than the 
Galactic field (Ambartsumian|1947] [1949). This definition, 
particularly the last part, distinguished associations from 
the other types of stellar group known at the time, open and 
globular clusters. Similar definitions were employed in the 


landmark studies by (1964), (1994), 
[Zeeuw et al.|( 1999) and|Lada and Lada|(2003). 


With the discovery of very young embedded clusters 
some studies argued that most embedded clusters and young 
star forming regions evolve into associations (e.g., 


2003; Zinnecker and Yorke|2007). In this context 
Lada and Lada|(2003) separated star clusters into embed- 


ded and open clusters depending on their association with 
interstellar matter. They argued that as embedded clusters 
emerge from their natal cloud they either survive as gravita- 
tionally bound open clusters or expand as unbound associ- 
ations. From this a more physical definition of associations 
emerged, one in which associations are specifically gravi- 


tationally unbound (e.g., [Gouliermis|2018 2020), 


a definition in harmony with that of a star cluster as being 


gravitationally bound (e.g., |Portegies Zwart et al.|2010). 


The difficulty with this physical definition is that, even 
for relatively isolated clusters and associations, it can be 
difficult to reliably determine whether a group of stars is 
gravitationally bound. This is not usually a problem for 
most classical OB associations since their stellar density 
is sufficiently low that there is no doubt that they are un- 
bound. However, it can be an issue for particularly young 
systems or where interstellar matter is present, since the 
precise mass and spatial distribution of this matter is much 
harder to determine than it is for stars. Furthermore, recent 
studies suggest star formation occurs over a continuum of 
densities rather than at either low or high densities (e.g., 
and 
therefore applying a density threshold to very young (still 
embedded) systems would be meaningless. 

For this reason we choose to define associations as 
groups of young stars with a stellar density lower than that 
of the Galactic field and that are not strongly associated 
with interstellar matter. The density definition separates 
associations from open clusters, while the interstellar mat- 
ter definition separates them from star forming regions and 
embedded clusters. This definition is advantageous because 
it is both purely observational and preserves the historical 
definition of associations. 

One note to this definition is that since associations are 
large and highly substructured they may (and often do) con- 
tain open or embedded clusters and star forming regions 
within their borders. The presence of star clusters as the 
nuclei of associations has been known since their original 
definition (Ambartsumian]|1949), while the the existence of 
star forming regions within their borders is evident from 
many of the nearest associations such as Sco-Cen or Orion 
OB1. These denser or younger regions are clearly part of 
a larger association, though they would not be classified as 
associations in their own right. 


1.3. This and past reviews 


This review summarises recent results on associations 
and our current knowledge of their properties and origins. 
It builds on a strong history of review articles on this topic 
dating from both the pre-Hipparcos period 
[1991), the Hipparcos era (during which many great discov- 
eries concering associations were made, e.g., 
1999) and the early Gaia era (Wright|2020). The wealth of 
studies concerning associations that have emerged over the 
last few years thanks to Gaia data have accelerated changes 
in our understanding of these systems, making this review 
very timely. 


2. Identifying members of associations 


Historically, the main obstacle to studying associations 
was the difficulty reliably identifying their members, partic- 
ularly the more-numerous low-mass stars that cannot easily 
be distinguished from older field stars. In this section we 
will discuss how the young low- and high-mass members 
of associations can be identified, the limitations of these 


techniques, the methods used to verify youth using spec- 
troscopy, and how samples of young stars are commonly 
divided into different associations or subgroups. 


2.1. 


There are numerous methods that can be used to iden- 
tify candidate young stars, depending on their effective tem- 
perature and age. Each of these methods requires different 
data, is effective under different circumstances, and has cer- 
tain biases that can limit the effectiveness of the method. 

Massive members of associations can be identified us- 
ing photometry. Following the reddening-invariant colour 


method pioneered by Johnson and Morgan|(1953), which 


maps out where stars lie in a colour-colour diagram depend- 
ing on their reddening, selected 
massive members of the cluster Westerlund 2 and surround- 
ing associations using the (u — g, g — r) colour-colour di- 
agram. Such diagrams show that reddened OB stars are 
located above the main stellar locus, and thus can easily 
be separated from the older field star population. Follow- 
up spectroscopy has demonstrated that this method pro- 
duces levels of contamination from late-type stars as low 
as 3% (Mohr-Smith et al.|2017). In regions of high extinc- 
tion (where u-band photometry is unavailable) the optical 
identification of OB stars becomes difficult. It is there- 
fore necessary to employ either near-IR or a combination 
of optical and near-IR photometry. 
used JH K, 2MASS photometry to select massive stars in 
Cygnus OB2, while (2018), 


and Quintana and Wright\(2021) used a combination of op- 
tical (Gaia) and IR (2MASS) photometry to identify mas- 


sive stars in the Galaxy within 4-5 kpc from the Sun. These 
samples are mainly contaminated by evolved massive or in- 
termediate mass stars (e.g. yellow and blue super-giants), 
but they are substantially free from RGB and AGB stars 
(see Section 2.3 in Zari et al.|2021). Finally, late B-type 
stars can live for +100 Myr, and can therefore be much 
older than the typical ages of associations (see SectionB.4); 
to definitely assign such stars to single associations it is 
necessary to combine photometric identification with other 
methods, such as kinematics. 

Historically, low-mass stars with ages of a few Myr have 
been difficult to identify, however dusty young stellar ob- 
jects (such as protostars or pre-main sequence stars with 
disks) can be selected through near- and mid-IR photom- 
etry, allowing for their identification with targeted surveys 
of individual star forming regions with telescopes such as 


Spitzer or Herschel (Gutermuth et al.|2009 
2017; Winston et al.|2020), or through all-sky surveys such 
as AKARI or WISE (Toth et al.[2014] Marton et al. 2016). 


These methods are effective, but heavily biased towards 
very young stars (« 5 Myr) that retain their disks, while 
disk-free candidate YSOs identified using IR photometry 
appear to be mainly contaminants (Manara et al.|2018). 
Identifying older stars and the substantial fraction of young 
stars without disks requires alternative methods. 


Identifying young stars 
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Fig. 1.— Colour-absolute magnitude diagram for stars in 
Orion based on Gaia EDR3. The black solid line is a 20 
Myr solar metallicity PARSEC isochrone 
[2017). For stars with G — Grp > 1, where the isochrone 
separates from the main sequence, isochrones like this can 
be used to separate the young PMS population above the 
isochrone from older field stars below it. 


Most low-mass members of associations are in the pre- 
main sequence (PMS) phase of stellar evolution. PMS stars 
are more luminous than main sequence stars of the same 
Teg and therefore stand out, above the zero-age main se- 
quence, in optical colour-absolute magnitude diagrams (see 
Fig. m). This provides an effective and efficient technique 
to identify young stars that requires only photometry and 
astrometry. The method is particularly useful for identify- 
ing low-mass stars (as low-mass PMS stars are more over- 
luminous at a given age than more massive stars) and isn’t 
biased towards other properties of the star (e.g., the pres- 
ence of circumstellar material or magnetic activity), how- 
ever some care is required when using it. For example, the 
extinction correction required by this method could intro- 
duce selection biases to the sample (see [Zari et al.[2018). 
Furthermore, the samples selected may be contaminated by 
unresolved multiple stars or evolved field stars with un- 
certain parallaxes, especially for PMS stars at intermediate 
masses that are less offset from the main sequence than low- 
mass stars. Finally, as Gaia parallax uncertainties increase 
as a function of magnitude, faint sources at large distances 
can be excluded when applying cuts on the relative parallax 
error 0,,/@. Nevertheless, following Gaia DR2 this ap- 
proach has been used in numerous studies to identify PMS 


stars younger than 20 Myr (e.g., 
2019). Attempts at 
identifying PMS stars older than 20 Myr (and younger than 
50 Myr) using this method have been made by |Kerr et al.| 
and|McBride et al.|(2021), who also propose differ- 


ent methods to estimate and correct for extinction. 
In the absence of reliable parallaxes, single epoch opti- 


cal or near-IR photometry can be combined to identify can- 


didate M-type stars or employed to deter- 
mine the size and spatial distribution of the low-mass PMS 
population of an association (see e.g. 
Armstrong etal O18). 


Single epoch photometry alone cannot definitively identify 
an individual star as a PMS star, thus samples selected based 
on this technique can include both background giants and 
foreground main-sequence stars of any age. 

Optical photometric variability is one of the defining 


characteristics of pre-main-sequence stars (Joy|1945 
ibig|1962). (2005) used variability combined 


with position in colour-magnitude diagrams and follow-up 
spectroscopy to select the low-mass young population in 
the Orion OBI association. This method is potentially bi- 
ased towards PMS stars with only high amplitude photo- 
metric variations, nevertheless current and future Gaia re- 
leases will facilitate all-sky searches for variable PMS stars 
and the use of variability as an additional selection criterion. 

PMS stars are magnetically active, which makes them 
bright X-ray sources. Studies of low-mass X-ray emit- 
ting young stellar populations of associations were carried 
out using data from the Einstein observatory (Walter et al.| 
1994), and the ROSAT all sky survey 
Walter et al.|2000). More recently, the X-ray observato- 
ries Chandra and XMM-Newton have been very effective at 
identifying large populations of young, magnetically-active 
stars, thanks primarily to their arcsecond-scale point spread 
function and high sensitivity (e.g., the Massive Young star- 
forming complex Study in Infrared and X-rays, MYStIX, 
which studied young stars in 20 Galactic massive star- 
forming regions, |Feigelson et al.[2013). However, the small 
fields of view of these observatories have meant they have 
preferentially targeted compact star clusters rather than dif- 
fuse associations (with the exception of more distant as- 
sociations, e.g., 2010). All-sky data from 
eROSITA overcome this issue 
and will allow members of many nearby associations to 
be identified, as shown by the preliminary study presented 
in Schmitt et al.) (2021) (see $6.3). As for the other tech- 
niques presented in this Section, although X-ray selection 
may weed out a lot of old stars, there will be many false 
positives remaining. For instance, the spin-down timescales 
of low-mass stars are ~50 Myr for solar-type stars and 
longer at lower-masses, so contamination by young field 
stars is probable. In addition, tidally locked, close bina- 
ries can maintain their X-ray activity indefinitely and appear 
*young' in the HR diagram due to the combined luminosity 
of the two stars. 

In summary there are many different methods for iden- 
tifying young stars. The advantages and disadvantages of 
each method, as applied to associations, are a combination 
of: whether the method is sensitive in the typical 1-50 Myr 
interval we are interested in; whether they offer any age 
discrimination within this range; whether they can be ap- 
plied to isolated stars; and whether the observational inputs 
are available over wide fields of view. Unfortunately, while 


each of the methods discussed here allows identification of 
young stars, they can all introduce contamination by older, 
unrelated field stars. These methods are therefore each in- 
sufficient in isolation to confirm the youth of stars. 


2.2. Confirming the youth of stars 


The secure identification of association members is in- 
extricably linked to confirmation of their age, with the most 
decisive indicators of youth coming from spectroscopy. Es- 
timates of spectral type or temperature can improve redden- 
ing determinations, allowing more accurate placement in 
the HR diagram (and this is often all that is required to con- 
firm the status of young high-mass stars), but for the more 
numerous lower mass stars there are more direct methods 
of determining youth. 

Lithium is an ephemeral element in the photospheres of 
low-mass stars; it is consumed by nuclear reactions on the 
PMS. Whilst a comprehensive understanding of the Li de- 
pletion process is still lacking, it is clear that Li abundance 
can be used as an empirical indicator of youth in a mass- 
dependent way (see Figure 2p. At low masses (< 0.7Mo) 
the presence of undepleted Li is a sure sign that a star is 
younger than about 20 Myr. Conversely, stars which ex- 
hibit total Li depletion must be much older, whilst interme- 
diate cases offer some age discrimination (effries|2014). 
Li is less diagnostic at higher masses because there is little 
PMS depletion and depletion on the main sequence is much 
slower. Using Li abundance as a tracer of young stellar pop- 
ulations in the nearest associations has a long history (e.g., 
Preibisch et al.|2002). Improvements in instrumentation 
are now allowing an extension of these techniques to wider 
fields and greater numbers of stars (e.g., Jeffries et al. 
quirement for reasonably high resolving power and signal- 
to-noise still limits this technique to associations closer than 
~ 1 kpc. 

Spectroscopic indications of gravity can be used to con- 
firm stellar youth at low masses and can resolve the con- 
fusion in absolute colour-magnitude diagrams caused by 
variability, reddening and binarity. PMS stars have lower 
gravities than MS stars of similar Tp, but for a given PMS 
age the difference is wider at lower masses, so these tech- 
niques are most sensitive for very young stars and/or PMS 


stars of low-mass (< 1M») (see Figure 8 of|López-Valdivia 
[et _al.|[2021 for an example). Surface gravity is difficult 


to measure directly, although its relative magnitude can be 
inferred from gravity-sensitive spectral indices in the opti- 
cal/near IR, such as atomic alkali, CaH or TiO molecular 
lines (Wilking et al./2005), empirically constructed spectral 
indices (Damiani et al.|2014) or the shape of the H-band 
peak (Scholz et al.|2009). Directly-measured surface grav- 
ities can be used to calibrate pre-MS evolutionary models 
(Olney et al.|2020) and directly infer the age of a low-mass 
star. In comparison to photometrically derived ages, this 
method does not underestimate the ages of stars on the bi- 
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Fig. 2.— EW of the Li 1 6708À absorption line for stars to- 
wards the y Vel cluster and Vela OB2 association. The bold 
diamonds show objects selected as members of the clus- 
ter/association on the basis of their Li abundance and posi- 
tion in the V vs V — I CMD (seeUeffries et al.[2014), while 
the open circles and open triangles (indicating upper limits) 
are sources not selected as members. The solid lines show 
theoretical predictions of the strength of this line at ages of 
10 (upper) and 20 (lower) Myr (Baraffe et al.|1998), while 
the dashed line indicates an empirical upper limit to the line 
strength at 50 Myr judged from observations of stars in the 
IC 2391/2602 clusters. The age discrimination of this Li 
feature is quite dependent on colour/mass; in particular, it 
is less sensitive at bluer colours/higher masses and may not 
offer a complete selection of members where the associa- 
tion stars are old enough to have depleted their Li in certain 
mass ranges (e.g. in the mid M-dwarfs in this example). 


nary sequence. 

Young stars with a protoplanetary disk are likely to be 
accreting gas. Classical T Tauri stars can be easily identified 
through having strong and broad Ha emission, even in low- 
resolution spectra (White and Basri|2003). These strong 
signatures are limited to just a fraction of the low-mass pop- 
ulation of an association, reducing from ~ 60% at a 1-2 
Myr to < 5% at 10 Myr (Fedele et al./2010). The comple- 
mentary population of “weak-lined” T-Tauri stars can still 
be identified from their weak, chromospheric Ha emission 
(along with several other lines, such as Ca II H & K). How- 
ever, these relatively weak activity indicators can persist for 
a few hundred Myr in low-mass stars so such samples in- 
evitably suffer contamination from field stars and active bi- 
naries, in a similar way to X-ray-selected stars. 


2.3. Assigning young stars into groups 


When a sample of candidate young stars has been iden- 
tified they can be divided into groups based on their distri- 
bution in position (plane of the sky and distance) and kine- 
matics (proper motion and radial velocity). This may be 
performed after the youth of these stars has been verified 


(see Section 2.2) or to confirm youth via the corroboration 
that a group of stars is spatially or kinematically coherent 
(the clustering of stars, both spatially and kinematically, de- 
creases as stars age, and therefore can be used as a youth 
proxy). A combination of kinematics and spectroscopic in- 
dicators of youth is increasingly used to get a sample of 
the low-mass population with minimal contamination (e.g., 
Wright er al 019). 

For compact clusters, stars can be grouped together us- 
ing just plane-of-the-sky information, often starting with 
a curated membership list (that removes most of the field 
contaminants) and then identifying the group using tech- 
niques such as the minimum spanning tree (Gutermuth et al.| 


2009), stellar density maps (Carpenter 2000||Megeath et al. 
2016), various mixture models (Kuhn et al.|2014), or us- 


ing other plane-of-the-sky clustering techniques 
al.|2019). Some of these approaches require a-priori as- 


sumptions regarding the number of different substructures 
present in a given population, and therefore the results of 
such methods should not be used to study the structure of 
such regions without understanding the biases at work. 

For low-density structures such as associations, assign- 
ing stars into groups cannot be done using just sky posi- 
tions. Since associations have relatively small internal ve- 
locity dispersions they continue to form coherent structures 
in velocity space even as they disperse. In the Hippar- 
cos era the identification of massive members of associa- 
tions relied strongly on proper motions (e.g.|de Zeeuw et al.| 
1999), sometimes combined with positions and parallaxes 
(de Bruijne [1999] [Hoogerwerf and Aguilar 1999), or ra- 


dial velocities (Rizzuto et al.|2011). Multi-epoch spectra 
are sometimes required to reliably identify binaries (that 


can introduce discrepant RVs due to the orbital motions), 
but even without multiple epochs, the broadening of an RV 
distribution by binaries can be estimated and assigning high 
membership probabilities based on a single-epoch of spec- 
tra is possible for most single and unresolved binary stars 
(es. González and Alfaro[pOT7), es- 
pecially in conjunction with other spectroscopic tracers of 
youth (see Section 2.2). 

With distances becoming increasingly common from 
Gaia, studies of the clustering and grouping of stars are now 
performed similarly in a higher-dimensional space (com- 
bining position on the sky, parallax, proper motion, and ra- 
dial velocity when available). As these multi-dimensional 
data combine different quantities (position and velocity) 
with different units, some prescription needs to be consid- 
ered to scale the data to give each dimension an appropriate 
weight. With higher dimensionality it is possible to select 
association members without initial identification of candi- 
date young stars (as is commonly needed for 2D cluster- 
ing). Several clustering algorithms have been used for such 


a purpose, including DBSCAN (Zari et al.12019| 
2021), and UPMASK (Canrar-Caudin et al 20195), Gaus- 
sian mixture models (Kuhn et al.|2020), as well as other 


custom codes. The typical inputs for these are the minimum 


number of stars in a group, the typical group dimensions, 
and the density threshold required to detect groups (though 
some codes are able to automatically adjust this threshold 
to detect both compact and extended populations). 

In recent years, with a large number of clustering codes 
and approaches available, as well as the variance that would 
result in using the same clustering algorithm tuned with 
different critical parameters, or applied to a differently- 
selected initial source list, it is common to see different 
groups identified within the same region of space in differ- 
ent studies (e.g.,[ounkel ef al 2018 
[et al.|2021). With the greater degree of attention given to 


extended associations (that are not as simple to define as 
compact star clusters) some confusion may arise regarding 
which of these overlapping groups should be considered as 
fundamental. As such, it is important to keep in mind that 
these algorithms only trace the underlying density distribu- 
tion of a given population with a greater or lesser sensitivity 
to the structures at a particular scale, often requiring a trade- 
off between the ability to identify finer substructure vs the 
ability to study the population as a whole. 


3. Properties of OB associations 


Associations come in diverse shapes and sizes and can 
appear, at first sight, a rather heterogeneous group. In this 
section we review the general properties of associations, 
with reference to notable individual examples, and the sim- 
ilarities and patterns that exist. 


3.1. 


The sizes of associations range from tens to several hun- 
dreds of parsecs (e.g., 
can vary between stud- 
ies due to different methods for defining the borders, differ- 
ent membership of the systems included, and on which sys- 
tems were included within the sample. For example, 


many and Stencel| (1992) measure a mean OB association 
size of 137 + 83 pc, while |Mel’Nik and Efremov) (1995) 


measure an average diameter of ~40 pc. Studies of extra- 
galactic OB association populations reveal broadly similar 
size distributions to Galactic associations, depending on the 
resolution of the data used. For example, [Lucke and Hodge| 
list 122 OB associations in the Large Magellanic 
Cloud (LMC) with sizes of 15—150 pc and a mean of 80 pc, 


while|Gouliermis| (2018) used the better-resolved catalogue 
of associations in the LMC from|Bica et al.|(2008) to mea- 
sure an average size of 30 pc. 

Associations increase in size with age (e.g., 
[1964), as would be expected for gravitationally unbound 
and expanding systems. Furthermore, associations con- 
nected to nebular material are typically smaller 
mis|2018). If the association with nebular material is an in- 
dication of relative youth, then this also supports a picture 
whereby associations increase in size as they age. 

Associations exhibit considerable internal substructure. 


Size and internal structure 


This includes the presence of subgroups with different ages 


and kinematics (Blaauw|1964||Garmany and Stencel 
young or open clusters (Ambartsumian||1949), or bright 


central concentrations (Ivanov) 1987||Mel’Nik and Efremov 
1995). In some cases this substructure is very clear, while 


in other systems it can be revealed using structural diag- 


nostics that quantify physical substructure (e.g., 
and Whitworth|2004| Wright et al.|2014). 


The most well-studied associations, such as Sco-Cen or 
Orion OBI, have been sub-divided into “OB subgroups” 
historically based on their on-sky distribution (e.g..[Blaauw| 
divided Sco-Cen into the subgroups of Upper Sco, 
Upper Centaurus-Lupus and Lower Centaurus-Crux). This 
physical substructure is often correlated with kinematic 


(Wright et al.\|2016}|Wright and Mamajek|2018) or tem- 
poral substructure (Pecaut and Mamajek;2016) that suggest 


these subgroups are real and not chance over-densities. 
Recent studies have combined on-sky positions with par- 
allax, proper motion or radial velocity information to sub- 
divide associations based on their 5- or 6D spatial and kine- 
matic structure (e.g., 
20190} [Berlanas et aL.|2019). For Da 
) and (2021) studied the Sco-Cen 


association, finding a wealth of sub-groups with different 
spatial and kinematic properties, which do not adhere well 
with the classical sub-division of the association. 
studied the young stars in the Vela 
OB2 association and found that the region is highly sub- 
structured, with seven main kinematic groups overlapping 
on the sky and whose formation is spread over ~ 35 Myr. 

These physical subgroups often show distinct age differ- 
ences, supporting the correlation between spatial, kinematic 
and temporal substructure (e.g., Zari et al.[2019). Further- 
more, they often extend beyond the bounds of the classi- 
cally defined OB associations (usually based on the distri- 
bution of OB stars) to expose older parts of the association 
(e.g. |Canrar-Gaudin et al [2019b). 

Many associations also contain open or embedded clus- 
ters within their borders that have been considered part of 
the association. In many cases their ages and kinematics 
are consistent with being related to the association, such as 


the y Vel cluster in Vela OB2 (Jeffries et al.[2014), p Oph 
in Sco-Cen 1991), or NGC 2353 in CMa OBI 
(Fitzgerald et al.|1990). In some cases however the open 


cluster has been shown to either have a very different age or 
significantly different kinematics, or both, suggesting that 
the cluster did not form as part of the association and is now 
just projected against the association (e.g., the NGC 2547 
cluster projected against Vela OB2, |Sacco et al.|2015). 


3.2. Kinematics, velocity dispersion and virial state 


Kinematic studies of associations were historically lim- 
ited to their most luminous members and used either Hip- 
parcos proper motions or radial velocities obtained from in- 
dividual spectra. The former were of limited precision for 
detailed kinematic studies (allowing membership to be con- 


strained, but rarely probing internal kinematics), while the 
latter were scarce and strongly affected by unresolved bina- 
rity. It wasn’t until the availability of multi-object spec- 


troscopy (e.g., [Preibisch et al.|D002) |Fitrész et al.|2008 
2016) and, later, Gaia proper motions (e.g., 
Wright and Mamajek|2018| |Kounkel et al.|2018) that de- 


tailed kinematic studies of the (more numerous) low-mass 
members of associations became possible. 

From kinematic studies using either radial velocities, 
proper motions or both, estimates of the velocity dispersion 
of associations can be calculated. Combining the velocity 
dispersion with estimates of the total mass of the associa- 
tion, the virial state can be derived. 

The accurate calculation of the velocity dispersion from 
measurement of the radial velocity or proper motion for in- 
dividual stars requires that the measurement uncertainty and 
the impact of unresolved binaries be accounted for. The for- 
mer can be easily quantified, but requires some modelling 
as the uncertainties are often significantly heteroskedastic 
or do not follow a normal distribution (e.g., 
[2015]find that the uncertainty distribution of Gaia-ESO Sur- 
vey radial velocities better follows a Student's t-distribution 
with extended tails). Unresolved binaries can significantly 
affect the instantaneous measure of velocity from single- 
epoch radial velocity surveys, particularly for high-mass 
stars that are EN found in binary or multiple systems 


(e.g., D TIS and require full modelling (e.g., 


The impact of unresolved binarity 
on EE proper in has not been fully explored, 


and while smaller than the effect on radial velocities, is still 
potentially significant (Jackson et al.|[2020). Proper mo- 
tions can also become unusable at large distances due to 
the distance-dependence of velocities derived from them, a 
problem that does not affect radial velocities. 

The corrected velocity dispersions of associations can 
show considerable variation, but are typically a few kilo- 
metres per second along each axis. For example, Ward and| 
studied 18 associations using Gaia DRI 
proper manone and measured velocity dispersions of 3— 
13 km s^! with a median of 7 km s^!, and 
[Dambis| (2020) measured an average velocity dispersion of 
4.5 km s^ ^ from 28 associations studied with Gaia DR2. 
Small differences between studies are predominantly due 
to differences in the association membership lists, many of 


which date from the 1980s (e.g., Humphreys|1978 
and Humphreys||1989) and need revisiting and updating 
with modern data (e.g., Quintana and Wright|202 


Immediately apparent from the earliest kinematic stud- 
ies of associations was that their velocity distributions were 
far from Gaussian and correlated with spatial position (e.g., 
[Tobin eraL]D009). Later studies us- 
ing proper motions (that probed more than one kinematic 
dimension) found that the velocity distributions of associa- 
tions were also highly anisotropic, with velocity dispersion 
ratios between the two proper motion axes of at least ~1.5 


(Wright et al.\2016}|Mel’nik and Dambis|2017 


Kruijssen|2018) and up to ~6 (Melnik and Dambis|2017). 


3D kinematic studies of associations using proper motions 
and radial velocities are rare, but similar levels of anisotropy 
have been found (e.g., Wright and Mamajek|2018). 

The anisotropy and non-Gaussianity observed within as- 
sociations has been attributed to kinematic substructure, 
groups of stars in the same area of space with similar kine- 
matics, but within a wider group with more diverse kine- 
matics. An example of this in Cyg OB2 from 
is shown in Figure B| where the proper motion vec- 
tors are coloured according to their position angle on the 
sky to highlight the kinematic substructure. This kinematic 
substructure has since been observed in Wight ad Mamaiel 2018 Kou mE 
et al.[2018] [Cantat-Gaudin et al. MOM et al.|2019b). "iere deren E 
many efforts to quantify the kinematic substructure. For 
example, Arnold et al.| (2020) analysed the kinematic sub- 
structure in Cyg OB2 and found a strong spatial/velocity 
correlation on sub-pc scales, but no correlations or struc- 
tures were found on larger scales. 

This kinematic substructure represents subgroups of 
stars within the association, and using sufficiently-accurate 
kinematic data the subgroups can be separated. There 
is growing evidence that when association subgroups are 
identified using kinematics that the velocity dispersion 
anisotropy of the subgroups is less pronounced (e.g., 
find velocity dispersion ra- 
tios of 1.0-1.6 for the subgroups in Vela OB2), which may 
indicate that the true association subgroups have isotropic 
velocity dispersions or might highlight a bias introduced by 
using kinematics to identify subgroups. 

The virial state of the association can be estimated by 
combining the measured velocity dispersion with estimates 
of the total stellar mass of the association (often extrapo- 
lated from a subset of the total population and therefore re- 
quiring assumptions about the form of the mass function), 
though see (2016) for caveats. Given 
that associations have a low stellar density compared to 
gravitationally bound open clusters, it is unsurprising that 
all associations to date have been found to be super-virial. 
The ratio between the virial mass and the stellar mass was 
found by [Mel nik and Dambis| (2017) to range from 10 to 
1000, with a median of ~50. Assuming there are no ad- 
ditional forces acting on the members of the association as 
it expands, then the virial mass will increase as the associ- 
ation expands (since the virial mass is proportional to ra- 
dius, [Portegies Zwart and Jílková|2015). If the expansion 
is assumed to be linear then the virial mass for an associa- 
tion will increase linearly with time, and as the stellar mass 
remains approximately constant, the virial to stellar mass 
ratio will also increase with time. 


3.3. Expansion 


Given their low space density, associations should be 
unbound and in a state of expansion and may have been 


more compact in the past (e.g., 1947| |1949 


1952). This expansion can be inferred from the 
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Fig. 3.— Proper motions for 798 X-ray and spectroscopically selected stars towards Cygnus OB2 (including 16 O-type 
stars), with kinematic outliers removed. The vectors are coloured based on their direction of motion to highlight kinematic 
substructure. The grey box shows the border of the X-ray observations used to identify members and the colour wheel in 


the top-right corner shows the relationship between colour and position angle. Figure from|Wright et al.|(2016). 


observed correlation between radius and density (Pfalzner| 
[2009). However, the internal velocity dispersions for many 
associations are too small to explain their present-day size 
by expansion from a significantly more compact state (e.g., 
[Torres er al 2008), raising 
questions over this interpretation. 

Measuring the expansion of associations is, in principle, 
very simple. put forward the linear expan- 
sion model, which assumes that associations had an initially 
compact configuration and have expanded linearly (in time) 
from then. For nearby associations it is necessary to ac- 
count for virtual expansion caused by the radial motion of 
the association towards (or away from) the observer, which 
can cause a false expansion (or contraction) on the sky, even 
when no physical expansion exists (Blaauw|1964). This can 
be corrected for using radial velocities, either for the bulk 


motion of the association or for individ- 
ual stars (Wright and Mamajek[2018) 

Early kinematic studies struggled to find evidence of ex- 
pansion in associations. For example, both [Wright et al.| 
and {Arnold et al.| (2020) could find no evidence for 


expansion in Cyg OB2 from ground-based proper motions, 


while|Ward and Kruijssen\(2018) searched for expansion in 


18 OB associations using Gaia DR1 data, but could not find 
evidence for expansion in any of their targets. The availabil- 
ity of Gaia DR2 astrometry lead to expansion being mea- 
sured in some, but not all, systems. Using Gaia DR2 data, 
studied 28 OB associations and 
found evidence for expansion within 6 of them. [Ward et al.| 
studied the kinematics of 110 OB associations using 
Gaia DR2 and argued that their properties were not consis- 
tent with expansion from a single, compact configuration, 
but were more consistent with originating from a highly 
substructured velocity field. 

More recent studies have found that if associations are 
divided into subgroups based on their spatial and kinematic 
substructure then these subgroups often show evidence for 
expansion, even when the whole system does not. For ex- 
ample, divided the Orion OBI asso- 
ciation into subgroups using spatial and kinematic informa- 
tion and found evidence for the expansion of the Orion D 


subgroup. |Cantat-Gaudin et al.|(2019b) divided the young 


stars of the Vela OB2 association into 7 subgroups and 


found all of them to be expanding. (2020) 


conducted a spectroscopic study of the centre of the Vela 
OB2 association to obtain radial velocities and found that, 
once the members of the Gamma Vel cluster were removed 
from the sample, the Vela OB2 association was expanding 
in all three dimensions (see Figure f}. 

For associations where expansion has been observed, 
the expansion is usually anisotropic, even when subgroups 
are identified using kinematics. Of the three OB associa- 


tions found to be expanding by|Mel’nik and Dambis\(2017), 


two are significantly anisotropic, with only Car OB1 being 
consistent with isotropic expansion. 
found that all three subgroups of Sco-Cen exhibit 
strongly anisotropic expansion, with all expanding prefer- 
entially along the Galactic Y axis but not along the other 


two axes. |Cantat-Gaudin et al.| (2019b) and 
(2020) also both found that the expansion of the Vela 


OB2 association subgroups was strongly anisotropic. 

In younger systems such as star forming regions or em- 
bedded clusters (which may represent the precursors of ex- 
panded associations), there is also growing evidence for ex- 
pansion. found evidence for expansion 
in >75% of their sample of young clusters, including the 
Orion Nebula Cluster, which |Da Rio et al.) (2017) had also 
found was expanding (though had not). 
observed a clear radial expansion pat- 


tern in the A Ori cluster, while |Wright et al.| (2019) found 
strong evidence for expansion in NGC 6530, though like the 


older systems the expansion of NGC 6530 is highly asym- 
metric, with almost all the expansion occurring in the dec- 
lination direction. 

In summary, there is growing evidence that many, if not 
the majority of, associations have substructures within them 
that are expanding. Older studies were unable to identify 
this expansion due to a combination of inferior kinematic 
data, poorly-defined association membership lists, or be- 
cause they were searching for expansion across the entire 
system rather than within the subgroups. The expansion 
that has been measured is generally anisotropic, a trend 
which also extends towards younger and embedded clus- 
ters. 


3.4. Ages and age spreads 


Establishing the ages of association members and quan- 
tifying differences, gradients and spreads in those ages are 
a crucial part of understanding the formation and evolution 
of associations. The only star with a model-independent 
age is the Sun, and beyond that there is a hierarchy of age- 
determination methods, of decreasing accuracy (e.g., see 
[Soderblom|2010). For young stars, these range from model- 
dependent methods such as asteroseismology, the fitting of 
stellar evolutionary isochrones in the HR, colour-magnitude 
or the log g vs Tog (Kiel, or spectroscopic HR) diagrams 
and the depletion and diffusion of light elements, to more 
indirect or empirically calibrated methods such as using the 
time-dependence of accretion, rotation and magnetic activ- 


ity (Soderblom et al.|2014). 
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Fig. 4.— Cartesian position-velocity plots of stars in Vela 
OB2 with 3D kinematics from[Armstrong et al.|(2020). The 
best-fitting correlation gradients between position and ve- 
locity are listed and shown in each panel as solid lines, with 
dashed lines indicating the uncertainties. All three dimen- 
sions show a positive correlation between position and ve- 
locity, which is a strong indication of expansion. 


The ages of high-mass stars can be estimated from their 
positions in the HR or spectroscopic HR diagrams. Age 
discrimination is possible because O and early B-type stars 
evolve quickly. Accurate spectroscopic determination of 
the stellar parameters and a distance (for the HR diagram) 
are required. Typically these techniques have been de- 
ployed in more distant but rich associations, where the large 
number of O-stars can be used to make statistical infer- 
ences about ages and age spreads. For example|Wright et al.| 
(2015) used 169 OB stars in Cyg OB2 to infer an age spread 
from 1-7 Myr. This work was expanded by |Berlanas et al.| 
(2020), with the benefit of Gaia DR2 parallaxes, and evi- 
dence is presented for ‘bursts’ of star formation (or at least, 
bursts of O-star formation) at ages of 3 and 5 Myr. Such 
studies are difficult because observational errors in the stel- 
lar parameters lead to significant age uncertainties and the 
fidelity of these ages can be compromised by unresolved bi- 
narity and a genuine astrophysical spread in the HR diagram 
caused by rotation-dependent internal mixing and mass-loss 
rate uncertainties (for the highest mass evolved stars). The 
relative paucity of high-mass stars can be countered by ac- 
cess to the larger low-mass populations, which further allow 
a comparison of the ages of the high- and low-mass stars. 

The low-mass populations of young (< 10 Myr) clus- 
ters and associations ubiquitously show a large spread of 
luminosity around any mean isochrone in the HR diagram. 
The causes of this spread are long-debated (see 
2014), 
Observational uncertainties play a role (particularly deal- 
ing with reddening, accretion, variability and binarity), but 
cannot account for all or even most of the dispersion (e.g., 
[Reggiani et al.|2011). If the dispersion were attributed to 
age, it would equate to age spreads > 10 Myr in many clus- 
ters and associations, but it is equally possible that astro- 
physical scatter associated with varied accretion histories 
or differing levels of magnetic activity might play a signif- 
icant role Gully-Santiago et al.[2017). 
This dispersion is a poorly understood nuisance factor when 
searching for age gradients or differences between popu- 
lations within an association and hampers interpretations 
involving sequential, triggered or bursts of star formation 
(though see for an instance in the 
Orion Nebula cluster, where multiple star formation bursts 
separated by ~ 1 Myr may have been resolved). 

Despite these difficulties there has been a frenzy of activ- 
ity in the past few years, exploiting both wide-field spectro- 
scopic surveys and Gaia astrometry to probe the age struc- 
ture of local associations. In Orion, 
looked at the stellar populations along the Orion A molec- 
ular cloud, using near-IR gravity diagnostics to confirm a 
spread in radius and thus possibly some age spread, but 
found little evidence for any age gradient along the cloud. 


Gaia DR2 studies by |Kounkel et al.| (2018), 
(2019) and|Kos et al.|(2019) have dissected the wider Orion 


OB association using clustering algorithms in position and 
velocity phase-space (see Figure [5] for an example). The 
ages quoted for these sub-regions range from 1-21 Myr, 
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spread over dimensions of ~100 pc. That the temporal sub- 
structure exposed by such studies shows a correlation with 
the spatial and kinematic substructure in these associations 
provides strong evidence that this substructure is real and 
reflects the star formation history within the association. 

A similarly complex pattern is emerging in the Sco-Cen 
and Vela OB2 associations. In Sco-Cen the rich, low- 
mass population has been used to trace clear age differ- 
ences (on average) between the younger Upper Sco re- 
gion at ^» 10 Myr and the older Upper Centaurus-Lupus 
and Lower Centaurus-Crux regions at ~ 16 Myr 
land Mamajek D016). (2021) have used 
Gaia EDR3 astrometry and spectroscopic radial velocities 
to show that about half the Upper Sco population is in the 
form of a clustered population that existed in more com- 
pact configurations in the past, whilst the rest is more dif- 
fuse. They suggest star formation proceeded over about 10 
Myr in small groups that gradually dissolve; indeed the de- 
duced kinematic ages (see Section B.5) of the compact sub- 
groups are younger than their isochronal ages. In a similar 
exercise for Vela OB2, (2019b) found 
seven subgroups with ages from 8 to 50 Myr. Age was not 
strongly correlated with position, but was correlated with 
kinematics suggesting a more turbulent than sequential star 
forming history. See Section[4.2]for a more detailed discus- 
sion of the star formation history within associations. 

In all these works it is explicitly (or implicitly) assumed 
that the model-dependent ages are accurate. Whilst there 
can be some confidence in the relative ages (or at least 
the age order) of different groups of stars, their absolute 
ages are more uncertain. There is a long-standing problem 
(sce and ref- 
erences therein) that more massive stars tend to have older 
isochronal ages than their lower mass siblings by factors of 
~ 2. Whilst there are undoubtedly problems to solve in the 
high-mass stellar modelling, it seems likely, with the emer- 
gence of evidence from PMS low-mass eclipsing binary 
systems (Kraus et al.|2015| David et al.]2019) and lithium 
eie Deer oar pace, that fitting 
conventional low-mass isochrones underestimates the ages 
(and masses) of low-mass PMS stars by factors of 1.5-2 and 
that the adoption of models incorporating rotation, mag- 
netic activity and starspots may bring these ages into much 
closer agreement with those of high-mass stars (e.g., 
— high-mass stars can also 


appear younger in the HR diagram if they are born rotat- 
ing significantly rapidly (e.g., Ekstróm et al.|2012). These 
uncertainties must be considered when discussing proposed 
scenarios where low-mass star formation is affected or trig- 
gered by the birth or deaths of high-mass stars. 


3.5. Kinematic ages 


Kinematic ages are often cited as providing a model- 
independent age for a group of stars, since they do not rely 
on any stellar physics. They are calculated from the time 
the system needs to expand from an initially compact con- 
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Fig. 5.— Estimated stellar ages across the Orion OBI association from |Kounkel et al.| (2018). Left: ages derived us- 
ing spectroscopic effective temperature and bolometric luminosities in the HR diagram. Right: ages derived using just 
photometry in the colour magnitude diagram (distances assigned using the average distance to stars in each group). 


figuration at its current rate to reach its present size (e.g., 
Makarov|2007). Their validity as 
a model-independent age is based on the assumption that 
the group of stars has expanded, unhindered, from an ini- 
tially compact state to its current configuration, and that the 
expansion began at, or shortly after, the birth of the stars. 
Kinematic ages have been calculated for many dif- 
ferent associations and moving groups (e.g., 
[2014). Historically, these have often disagreed 
with isochronal ages calculated for those groups (e.g., 
Soderblom\2010). However, the combination of improved 
Gaia astrometry, a revision of the pre-MS evolutionary age 


scale (e.g., Bell et al.|2013), and improved kinematic age 
calculation methods (e.g., 2018} |Crundall 
[et a1. |2019) have resolved many discrepancies. For exam- 


ple, the 8 Pic moving group has recently been calculated 


to have a kinematic age of 18.3*13 Myr (Crundall et al. 
2019), in excellent agreement with its lithium depletion 


boundary age of 21 + 4 Myrs (Binks and Jeffries|2014). 
The reliability of kinematic ages is dependent on a num- 
ber of assumptions (see for a general 
critique of this method). Foremost amongst these is the as- 
sumption that the time at which the stars were closest to 
each other was when they formed. An example relevant to 
this is the A Ori system, whose isochronal age suggests it is 


much younger than its kinematic age (Kounkel et al.|2018). 
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If the isochronal age is correct, and not underestimated (see 
Section [3.4]for a discussion on the reliability of isochronal 
ages), it would imply that the stars in the A Ori group may 
have formed from material that was already expanding. If 
this is a common occurrence it raises serious questions over 
the validity of kinematic ages. Related to this is the as- 
sumption that associations formed in a compact configu- 
ration (see Section B.3] and in particular the evidence for 
anisotropic expansion), with equally serious implications. 
performed N-body simulations to 
test the validity of kinematic ages and how astrometric un- 
certainties and the different methods employed to measure 
kinematic ages can affect the results. Their results sug- 
gest that the traceback method will underestimate the age of 
the association (and overestimate its initial size), with ages 
converging to ~4 Myr, while the method of comparing ve- 
locity with position in a given dimension can also lead to 
significant uncertainties. Approaches such as the forward- 
modelling method used by [Crundall et al.) (2019) have the 
potential to address some of these issues, but could be com- 
putationally time-consuming for large systems. 


3.6. Distribution within the Galaxy 


Since the seminal work by (1953), OB 


stars and OB associations have been used to trace the spatial 
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Fig. 6.— Left: Density distribution of OBA stars from [Zari et al.|(2021) in the Galactic plane. The density is displayed in 
arbitrary units. The Sun is at X, Y = 0,0. The x-axis is directed towards the Galactic centre, and the y-axis in the direction 
of Galactic rotation. The z-axis is perpendicular to the plane. The dashed circles have radii from 1 to 6 kpc. The black 
crosses indicate the position of the OB associations listed in [Wright] (2020). The black lines show the spiral arm model 
from|Reid et al.|(2019). Top Right: Zoom-in of the OBA star density map shown on the left. The white stars represent YSO 
groups identified by (2021b). The black thin line shows the location of the Sagittarius-Carina arm. Bottom 
right: Density distribution of dust from (2019). As in the left plot, black crosses indicate the position of 
OB associations listed in [Wright] . The thick gray lines represent the approximate location of the Radcliffe wave 
(left) and the Split (right). 


distribution of young stars in the Milky Way, and to probe — includes star forming regions in Orion, Taurus, Perseus and 
the spiral arms of the Galaxy. In recent years many studies | Cygnus. The Split is argued to be a 
have linked the spatial distribution of young stars, young long arm segment linking the Local and Sagittarius-Carina 
clusters, and associations in the solar neighbourhood (d < spiral arms, and includes the Sco-Cen association. Figure[6] 
500 — 1000 pc) into larger structures. (bottom right) shows a schematic representation of the Rad- 


Bouy and Alves| (2015) studied the spatial distribution — cliffe wave and the Split, plotted on top of the dust density 


of O and B-type stars within 500 pc of the Sun. They map from|Lallement et al.|(2019). 
suggested that the distribution of OB stars in the solar Thanks to data from Gaia it has been possible to char- 


neighbourhood is described by stream-like structures called — acterise the distribution of associations and young stars be- 
*blue-streams'. Such blue-streams are associated with the yond the solar neighbourhood. estimated 
three largest OB associations within 500 pc of the Sun: and studied the 3D distribution of OBA stars in the Milky 
Orion OBI, Vela OB2 and Sco-Cen. The work of|Bouy and| Way disk within 4-5 kpc of the Sun. Figure [6] (left) shows 
was based on data from the Hipparcos satel- the distribution of the OB associations from the list pre- 
lite, and motivated (2018) to perform a follow-up sented in Table 1 of and references therein) 
study using Gaia DR2. confirmed that projected on the Galactic plane, plotted on top of the den- 
the 3D structure of star forming regions in the solar neigh- sity distribution of the filtered sample of OBA stars from 
bourhood is complex and filamentary, although they did [Zari et al. (2021). In general, the distribution of OB associ- 
not find evidence for the presence of the ‘blue-streams’ hy- — ations and OBA stars trace similar structures. These struc- 
pothesised by [Bouy and Alves| (2015). Indeed, subsequent tures can be identified as part of the Sagittarius (or some- 
studies have found that numerous associations in the solar times Sagittarius-Carina) Arm towards the inner Galaxy, the 
neighbourhood can be linked to two gaseous structures: the Orion (or sometimes Cygnus-Orion) Spur approximately at 
Radcliffe wave and the Split. The Radcliffe wave the Sun's position, and the Perseus Arm towards the outer 
isaco- Galaxy. There are however a number of important dif- 
herent and narrow structure around 2.7 kpc long, whose 3D ferences. The distribution of OBA stars shows a strong 
shape is well described by a damped sinusoidal wave in the — over-density corresponding to the Scutum-Centaurus arm 
plane of the Milky Way. The Radcliffe wave corresponds — towards the inner Galaxy (at approximately X ~ 2 kpc). 
to the densest part of the Local Arm of the Milky Way and Only two associations seem to be loosely linked with this 
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arm suggesting our census of OB associations is incomplete 
at such distances. Many associations also do not seem to 
correspond to any significant over-densities in the OBA star 
distribution: such associations might have too low a den- 
sity to appear on Fig. [olor may just be over-densities of OB 
stars in the sky. Finally, the mean distances to numerous 
associations derived in the literature do not correspond to 
the distances of the over-densities in Fig. [e] although they 
are towards the same lines of sight. This calls again for a 
revision of the census, membership, extent and distances of 


associations (see e.g., Quintana and Wright\2021). 


Finally, young stellar objects (YSOs) can also be used to 
trace Galactic structure. presented a cat- 
alogue of candidate YSOs, and found groups of YSO can- 
didates associated with the Local, Sagittarius-Carina Arm, 
and Scutum-Centaurus Arms. used the 
same catalogue and focused on a linear feature between 
Galactic longitudes | ~ 4?—18.5? including the star form- 
ing regions M8, M16, M17 and M20. Figure [6| (top right) 
shows the YSO groups associated with this structure com- 
pared to the OBA star density map of|Zari et al.|(2021) (Fig. 
[e]left). The structure traced by the YSO groups corresponds 
to a prominent over-density that is also visible in the OBA 
star map and mapped out by the previously-known Sagit- 
tarius, Scutum and Serpens OB associations in this region 
from[Wrighr](2020). The structure does not seem to be iso- 
lated, but appears to be a feature of the Sagittarius-Carina 
arm, similar to those observed in external galaxies. 


4. Discussion 


Historically, star clusters have been much better stud- 
ied than associations due to the observational bias that 
they are much easier to study and categorise as they are 
spatially much more compact, and suffer far less from 
back/foreground contamination. Arguably, this has led to 
a bias towards clusters being considered much more im- 
portant than associations as sites of star formation. The 
frequent use of the phrase ‘most stars form in clusters’ is 
usually supported by a reference to[Lada and Lada|(2003). 
However, Lada & Lada define a ‘cluster’ as a group of 35 
or more physically related stars whose stellar mass density 
exceeds 1.0M pc^? — a broad definition which says noth- 
ing about boundedness. A more correct statement, and one 
more consistent with Lada & Lada’s view, would be ‘most 
stars form at significantly higher densities than the field’, 
followed by ‘and around 10% of stars remain in bound clus- 
ters at significantly higher densities for > 10 Myr’. 


4.1. 


There are many theories as to how GMCs form, and how 
they turn a globally supersonic medium into star-forming 
*units'. For a detailed overview of star formation, and, in 
particular, how GMCs convert gas into stars, we direct the 
reader to other chapters in this volume, such as Chapters 1, 
5 and 7. For the purpose of this chapter the important point 
is that when we observe stars with ages of a few Myr we 


Formation and origin 


observe them at both relatively low densities in associations 
and at high densities in bound clusters. A key question in 
star formation is: why? 

At one extreme is the view that all stars form at high 
densities in bound clusters, but that the vast majority of 
these clusters are very rapidly dispersed by gas expulsion. 
The essence of this idea is that an embedded star cluster in 
virial equilibrium will expel the residual gas left over from 
star formation thanks to feedback from massive stars on a 
timescale of a few Myr. Since the star formation efficiency 
is generally quite low and therefore the majority of the mass 
is still in the form of gas, expelling this gas will dramatically 
reduce the potential of the cluster causing it to unbind and 


quickly disperse (e.g. 1980} Kroupa|2001 
land Bastian|2006). Associations are then the expanding 


remnants of one (or more) originally dense star clusterd!] 

It has become clear that the view of associations as the 
remains of (a small number of) expanded clusters does not 
fit the observational data. Recent studies have revealed that 
associations are more spatially extended than previously 
thought, and that they present a high degree of substruc- 
ture in physical space, kinematics, and age (see Section 
3). This suggests that associations are globally dynamically 
very yeu] (Parker et al.\|2014), since any form of sub- 
structure is easier to erase than to form — any encounters 
within a group of stars will erase substructure within a few 
crossing times, which can be particularly rapid in bound 
stellar systems (e.g., 
worth\2004} Parker et al.|2014). This strongly suggests that 
associations form with a similar spatial configuration as we 
currently observe them: over large volumes, with consid- 
erable substructure, at low average densities (though some 
stars may form at higher densities in subclusters), and most 
likely globally unbound. 

The kinematic and physical substructure observed in as- 
sociations reflects that of GMCs. Observations of cores 
within GMCs show stars often form at relatively low den- 
sities with significant substructure. The detailed large-scale 
sub-mm maps from Herschel show cores (unsurprisingly) 
following the complex gas structures. Two particularly 
good examples are Aquila (Könyves et al.|2010| see their 
fig. 1), and Orion B their fig. 
5). Similar structures are seen in very young class I/II 
stars traced by e.g. Spitzer (Gutermuth et al.|2009). We 
would expect the distributions of cores and stars we see over 
scales of several pc in e.g. Aquila and Orion B to be only 
marginally bound or unbound (e.g. from Larson’s laws, 
Larson|1981). Therefore, it is difficult to imagine these re- 
gions evolving into anything other than associations (with 


1For much of the 1990s and 2000s this was one of the dominant pictures 
of star formation and an underlying assumption in much work on clusters 


e.g.|Goodwin and Bastian|(2006). 


?Note that dynamical age and chronological age can be very different. The 
former is a measure of age in terms of crossing times, which, in unbound 
systems, is effectively infinite (some relaxation can occur on small scales 
and within substructures, and unbound initially high density regions will 
have encounters early-on that can erase substructures, but the point stands). 
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maybe some small subclusters). 


4.2. Propagation of star formation in OB associations 


Whilst complex and difficult to measure, the ages of stars 
and age distributions within associations provide clues as to 
how star formation has proceeded in different regions. 

The classical model for the formation of OB sub-groups 
in associations was proposed by 
(1977). This model predicts that ionising radiation and 
winds from massive stars terminate the star formation pro- 
cess and drive shocks in other parts of the parental cloud. 
New generations of massive stars are born, and the process 
is repeated. As a consequence, low-mass stars should be 
older than OB stars, since only OB stars are formed by trig- 
gering, while low-mass stars should form spontaneously 
through the cloud. Another popular model is based on the 
mechanism of radiation driven implosion (Sandford et al.| 
Haworth and Harries|2012). According to this model, low- 
mass stars should be younger than the OB stars (which 
initiate their formation), and one may expect to see an age 


gradient in the low-mass population (Preibisch and Zin- 
2007). However, from an observational perspective, 


there is currently no convincing evidence for either of these 
mechanisms. 

Being the closest OB association to the Sun (d ~ 
140 pc), Sco-Cen has provided an ideal laboratory for test- 
ing theories of the propagation of star formation (Preibisch| 


1998| 2002; |Preibisch and Zinnecker|2007). Fol- 
lowing the pre-Gaia study by |Pecaut and Mamajek\(2016) 
above, (2018) 


and Gaia studies mentioned 
and [Kerr et al.| (2021) proposed more complex models for 
the star formation history of Sco-Cen, which take into ac- 
count the increasing layers of complexity emerging from 
the data. combined gas observations 
and hydrodynamical simulations to study the formation of 
the Sco-Cen super bubble, and suggested the following sce- 
nario for the evolution of the association. Dense gas was 
originally distributed in an elongated cloud, which occu- 
pied the current area of the association. The star-formation 
events in Upper Centaurus Lupus and Lower Centaurus 
Crux lead to super-bubbles that expanded, surrounding and 
compressing the parental molecular cloud, and triggering 
star formation in Upper Scorpius. This scenario predicts 
the formation of kinematically-coherent sub-groups within 
the association that move in different directions, similar to 


the observed kinematics in Sco-Cen (Wright and Mamajek 
2018). (2018) also predict that young groups 


could also be found in regions containing older stars, and 
that several young groups with similar ages might form over 
large scales. This is consistent with what is observed in 
Sco-Cen and other associations. The ages derived by [Kerr] 
however seem inconsistent with the model 
proposed by Krause eral QOIS), argue 


that Krause's 'surround and squash' scenario might act in 
certain regions of the association, while in others a simple 
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model of sequential star formation might have been enough 
for the propagation of star formation. It could therefore be 
that multiple methods of star formation propagation have 
acted within a single association. 

The shell-like distribution of stars and gas in the Vela 
OB2 complex suggests an episode of triggered star forma- 
tion that shaped the gas into a shell, causing it to compress 
and form stars (Cantat-Gaudin et al.|2019a). This would 
indicate that the expansion of the IRAS Vela shell preceded 
the formation of the stars in Vela OB2, with the expanding 
motion of the gas imprinted on the stars that formed. The 
mechanism responsible for shaping the distribution of stars 
and gas would have to be energetic enough to not only com- 
press gas into forming stars, but also induce the expansion 
of the IRAS Vela shell, which|Cantat-Gaudin et al.|(2019a) 
suggest was most likely a supernova. 

A similar conclusion was drawn in Orion by 
and (2021), who independently 
concluded that a major feedback event occurred in the re- 
gion ~ 6 Myr ago. Such an event shaped the spatial dis- 
tribution and kinematics of the gas and young stars that 
are observed today, and possibly triggered the formation 
of some of the youngest stellar groups in the region. Al- 
though [Rounkel (2020) and|Grojischedl er aT 2071) stud- 
ied the motions of different groups as a function of time, 
they did not assess the presence of age gradients at the time 
such groups were formed (instead assessing them at the cur- 
rent time). This may be important to better explain the for- 
mation history of associations and potentially to evaluate 
the importance of triggering. The star formation history of 
the Orion OB association is however not completely un- 
derstood. The age distribution of the groups identified for 
instance by [Kounkel et al.|(2018) and[Zari et al.|(2019) has 
not produced a clear picture of the progression of star for- 
mation across the association. 

It is interesting to compare the large age spreads and 
seemingly complex star formation histories of associations 
with those of clusters. In particular, the association Cyg 
OB2 and the young massive cluster Wd 1 are both a few 


x10* Mc with an age of ~ 5 Myr (Brandner et al.| 2008 
Wright et al.|2015). However, Cyg OB2 seems to have a 
significant age spread (Wright et al.|2015), whilst Wd 1 


seems to have formed in a single burst of star formation 
(Brandner et al.|2008). Presumably the very different star 
formation histories are telling us about the properties of the 
(presumably fairly similar mass) GMCs from which they 
formed. The reason why some GMCs produce massive OB 
associations, while other GMCs of similar mass produce 
massive clusters is, however, still far from clear. 


4.3. Dynamical evolution and dispersal 


OB associations are massive, extended stellar structures 
whose low density suggests they are globally unbound. 
Their expansion rates are sufficiently low (or are limited 
to one or two dimensions) that they remain co-moving for 
some length of time and form elongated structures. Many 


elongated, and relatively old, association-like structures 
have been discovered in recent years using Gaia data. 
reported the discovery of a 
= 17 Myr population, clustered in proper motions, but fila- 
mentary and extending for ~90 pc in physical space, which 
they interpret as a relic of star formation in a molecular 
cloud filament. An analogous structure was also reported in 


Vela by [Beccari et al.|(2020), extending ~260 pc and with 
an approximate age of 35 Myr. |Kounkel and Covey|(2019) 


also identify a number of filamentary structures orientated 
parallel to the Galactic Plane and with ages up to ~100- 
200 Myr (young enough that they have not completed a 
full Galactic orbit) — co-moving, filamentary and extended 
groups much older than this appear difficult to identify 


using currently-available data and methods (Kounkel and 
|Covey|2019). 


These groups are too young to have been significantly 
affected by tidal forces and therefore their extended struc- 
ture is likely to be mostly primordial. As associations are 
thought to form from filamentary molecular clouds, the re- 
sulting young population should retain a similar string-like 
morphology, which would explain the observation of these 
extended structures. Furthermore, since the dissolution of 
associations is not completely isotropic (see Section p.3). 
and in some cases shows preferential expansion in the X and 
Y directions of the Galactic Plane (e.g..[Wright and Mama-| 
[jek[2018), this may explain the origin of such extended and 
filamentary structures. 

As a population of stars is orbiting around the Galaxy, 
its velocity dispersion increases due to dynamical evolution, 
tidal forces (such as from passing stars, nearby GMCs, and 
the Galaxy itself) and shear due to Galactic rotation. Over 
time stars in the group will slowly drift away and eventually 
most stars in the group will have dispersed into the galactic 
field population. Some of the stars may still remain as a co- 
herent group for several tens if not hundreds of Myr, with 
the populations that were originally more massive manag- 
ing to retain a larger fraction of their members as co-moving 
for longer periods of time. iden- 
tified a large number of co-moving groups in the Milky Way 
using Gaia DR2 data. They found that the number of stars 
within a co-moving and coherent group decreases with time 
according to a power law. 


5. Implications 


Over the last decade there have been many changes to 
our view of OB associations and their origins, and these 
have significant implications for our understanding of the 
star and planet formation processes. 


5.1. At what densities does star formation take place? 


A key question for star formation studies in galaxies is in 
what sort of environments do stars typically form? Do most 
stars form at the fairly low stellar densities observed in as- 
sociations? Or do they often form at higher (cluster-like) 
densities and disperse? Or are both important (e.g. [Besser] 


15 


2010} 2012)? This is a question with poten- 


tially significant implications for star formation: do stellar 
systems form in cores that are essentially ‘isolated’, or do 
dynamics and encounters play an important role in the for- 
mation of stars and mass assembly? 

The distinction is typically made as to whether a partic- 
ular region is a ‘cluster’ or an ‘association’ from its total 
energy: clusters are bound, and associations are unbound 
(e.g. 2011). However, from the 
point of view of star formation, multiple systems or plan- 
ets, a more useful distinction is if stellar systems have ever 
‘known’ about other systems. This relies on their ‘den- 
sity history’ — in particular, the highest density environment 
they have ever been in, and for how long. 

There is a 'critical density in a star forming environment 
of ~100 Mo pc? above which star forming cores interact 
with each other. This comes roughly from how closely one 
can pack 0.1 pc radius cores with a lifetime of ~0.2 Myr and 
velocity dispersion of ~1 km s^! - 
above this density cores will interact with each-other while 
forming stars and not be ‘isolated’. Even if a region at this 
density is unbound, encounters will occur before it expands 
(unless the expansion velocity is 10s km s71). 

If, as seems clear from the data on associations, most of 
the stars in associations have always been at low density, 
then these systems will be ‘pristine’. That is, they have not 
been altered by external encounters/irradiation (however, 
secular effects may be important) and it will be important to 
compare the properties of such stars, together with their cir- 
cumstellar environments and planetary systems with those 
found in dense clusters. 


5.2. Do binary and multiple systems evolve differently 
in clusters and associations? 


A significant fraction of stars are thought to form in 
binary and higher-order multiple systems (see reviews by 
Kraus\2013). Binary and multiple systems can be changed 
or destroyed by secular decay (for higher-order multiple 
systems) and by external encounters. Secular decay should 
be environmentally-independent (although small external 
perturbations could play some role) and occurs on order of 
the multiple system's crossing time (typically 10s to 100s 
kyr). For example, a triple system would usually be ex- 
pected to decay to a (hardened) binary and a single star 


(e.g.|Goodwin and Kroupa|2005| Reipurth et al.|2010). The 


effect of close encounters depends on both the frequency 


and energy of encounters (Heggie|1975 1975 
1995) — both of which are density-dependent (though it is 


generally a very stochastic process, 
[2012). The more encounters and the more energetic they 
are, the more systems we would expect to be destroyed (or 
at least significantly altered). 

The effect of both secular decay and close encounters 
is therefore to decrease the multiplicity fraction (MF) of a 
group of stars over time. Unfortunately, these two effects 


are extremely difficult to disentangle as they both lead, in a 
very stocastic way, to a reduction in the MF. 

Observations seem to show that multiplicity decreases 
with age. The MF for Class 0 YSOs is close to unity 
et al.[2013] [Tobin et al.|2019), but by the time these sources 
reach the Class I and Class II stage the MF is closer to 
what is typically observed in the field (see 
Reipurth et al.[2014). However, the simplify- 
ing assumption of a single initial population of multiple sys- 
tems which are then altered by dynamics (e.g.[Kroupa|1995] 
seems difficult to justify a priori. This is because if 
the star formation environments are different enough to pro- 
duce an association versus a cluster, one might well expect 
that the initial multiplicity properties of stars would be dif- 
ferent. This leaves us with attempting to reverse-engineer 
from current populations the stochastic history of secular 
decays and external encounters on potentially different ini- 
tial multiple populations which could give rise to the cur- 
rently observed populations. 

From our earlier discussions that most stars in associa- 
tions seem always to have been at low density, while stars in 
clusters we know have spent a significant fraction (if not all) 
of their lives at high density, we would expect the multiplic- 
ity properties of clusters and associations to be different as 
the encounter/irradiation histories of stars within them were 
different. There is some evidence supporting this view. In 
the low-density Taurus Molecular clouds the multiplicity 
fraction is very high, ~66—75%, in excess of what is found 
in the field by a factor of ~2 for stars of a comparable mass. 
In contrast, in the dense Orion Nebula Cluster (ONC), the 
multiplicity fraction at separations of 67.5—675 AU is con- 
sistent with the field (Reipurth et al.|2007). 

The dependence of the MF on both local stellar density 
and binary/multiple system separation is also interesting, 
with wide (7 1000 AU) binaries particularly sensitive to en- 
vironment. For example in the low-density Taurus Molecu- 
lar cloud and Cygnus OB2 association there are many such 
wide binaries 
[2020), while in the dense ONC no binaries with separations 
over 1000 AU have been conclusively identified 


1999). These observations can help place limits on 
the density history of such binaries. For example, 
(2018) showed that the high fraction of wide binaries 


amongst massive stars in Cyg OB2 was inconsistent with 
them being born in a dense cluster with other massive stars, 
because wide binaries in clusters will be destroyed or hard- 
ened by the presence of another massive star. 

For close binaries (<100 AU), all nearby low-density re- 
gions that have been studied appear to show a MF that is 
very similar to the field for separations of 60-100 AU (with 
only a tentative hint of a density dependence) but an ex- 
cess of systems < 60 AU by a factor of roughly two 
[et al.|/2012). Separations of < 60 AU are interesting as 
these systems are ‘hard’ in all local environments and so 
should survive into the field. So an excess in local regions 
over the field of a factor of two suggests some environments 
must under-produce fairly close systems. The obvious sus- 
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pect for under-production are clusters. The ONC is the only 
dense cluster close enough for detailed multiplicity studies, 
and it has been difficult to resolve close (< 60 AU) com- 
panions in the ONC. A recent search for close companions 
in the ONC using high resolution imaging has shown a sim- 
ilar excess in the MF at 20—60 AU as in low density regions: 
twice the field (Duchéne et al.12018). However, this effect 


may be limited to particular primary mass ranges 
et al.|2019|2021). It is also problematic to extrapolate from 


a sample of one cluster. 

At intermediate separations (100-1000 AU) the picture 
is still unresolved. There is evidence that stars in the denser 
regions of the Orion Molecular Clouds have a higher MF at 
separations of 100-1000 AU than the regions with lower 
stellar density [robin er al 022). 
Thus, despite being presumably much more sensitive to 
destruction via encounters than low density populations, 
clustered regions can seemingly produce a larger fraction 
of companions at some separation ranges. This could be 
primordial, it could be due to the tightening of previously 
wider binaries, or it could be due to the loose capture of 
stars into multiple systems Moeckel and Clarke|(201 1). 

In summary it is currently impossible to say with any 
certainty if or how birth multiplicity varies with environ- 
ment (or if it is universal) and how it is processed differ- 
ently in different environments). In the future, in light of 
the larger census of YSOs that can now be assembled with 
a variety of different surveys such as Gaia, and with new 
facilities coming online that offer increasingly better spatial 
resolution, it may be possible to further improve multiplic- 
ity statistics in a larger number of star forming regions to 
better determine the role of the environment in the under- 
lying processes governing both the initial multiplicity of a 
population and its secular decay and the effect of external 
encounters. 


5.3. Do protoplanetary disks and planetary systems 
evolve differently in clusters and associations? 


There are two external effects that can alter the proper- 
ties of planetary systems: encounters and irradiation. En- 
counters can affect both disks and (N-body) planetary sys- 
tems, while irradiation only affects disks. The stellar den- 
sity of the environment surrounding a young star (and the 
evolution of that density) determines both of these effects 
through the likelihood of close encounters and the local ra- 
diation and particle field. 

In dense clusters, disks may be photoevaporated by EUV 
radiation from nearby O-type stars (typically when within 


1 pe, e.g. [Winrer er aL.]2018] [Cai er al. 019 
2021), while even a brief period (« 1 Myr) at high density 


is enough to alter planetary system properties (see 
2020| and references therein). At lower densities in asso- 
ciations these effects should be different and reduced in 


comparison (e.g. |Laughlin and Adams Clarke and 
1993| Adams et al.|2006| Wright et al.|2012. 
gies Zwart and Jílková|2015 2015 


et al.[2019} Winter et al.|2019b). 


Current observational evidence is limited, but suggests 
that disks are affected by their environment. |de Juan Ove- 
(2012) found evidence for a decrease in the proto- 
planetary disk radius at very high cluster densities, but their 
sample of sources at such high densities was small. 
(2018) made a similar discovery in the ONC, when 


compared to lower-density star forming regions. 
(2016) found that the fraction of stars with protoplan- 


etary disks is inversely proportional to the local ultraviolet 
flux in Cyg OB2, suggesting disk dissipation is dominated 
by photo-evaporation and therefore strongly dependent on 
environment. 

The extreme mass loss rate in the protoplanetary disks 
in the ONC from photoevaporation from 0! Ori C 
has raised the question of how these disks 
can survive for a prolonged period of time to be currently 
detected, without invoking implausibly large initial disk 
masses or unrealistic ages for 0! Ori C, given the typical 
ages of other stars in the vicinity. However, the age gradi- 
ent of stars in the ONC may provide 
a solution to the problem: older stars have had a greater 
shielding from photoevaporation through the surrounding 
dust, allowing them to circulate throughout the cluster and 
migrate away from 6! Ori C. On the other hand, the most 
strongly irradiated stars are the youngest, with the greatest 
mass reservoir remaining (Winter et al.|2019a). 

Further measurements of protoplanetary disk frequen- 
cies in OB associations will be needed to fully understand 
the influence of environment on disk evolution. To under- 
stand how the environment affects the properties (or fre- 
quency of) planetary systems will require large samples of 
planetary systems in star clusters and OB associations to 
be compiled. The current evidence and our understand of 
planet formation does suggest that the properties of plan- 
etary systems should depend on environment, though the 
exact form of this dependency could take many forms. 


6. Future Prospects 


Many upcoming facilities and instruments have the po- 
tential to significantly advance our understanding of associ- 
ations. Here we discuss these observatories, as well as the 
prospects of future Gaia data releases. 


6.1. 


The data from the Gaia satellite have allowed for signif- 
icant advancements in our understanding of associations. 
Thanks to Gaia's exquisite astrometric precision, it has 
been possible to identify low-mass members of nearby as- 
sociations (within 500 pc) and to study in detail their distri- 
bution in 3D physical space, their kinematics and ages. 

With Gaia EDR3 the precision of parallaxes and proper 
motions has improved on average by 2096 and by a fac- 
tor of 2 respectively, compared to DR2. Gaia DR3 (ex- 
pected mid-2022) will further provide new and improved 
radial velocities (for ~30 million stars out to G ~ 15 — 


Future Gaia data releases 
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15.5 mag, with precisions «1 km s^! for the brightest 
stars) and astrophysical parameters for sources based on 
the BP/RP prism and RVS spectra (for the full content 


of Gaia DR3 see https: //www.cosmos.esa.int/ 
web/gaia/release). Gaia DR4 will feature parallaxes 


that are more precise by a factor 1.7 with respect to Gaia 
DR2, while for proper motions the gain will be a factor of 
2 (Brownp021). 

The parallax precision improvement from future Gaia 
data releases, combined with a better understanding and 
treatment of the systematics, should allow the identification 
of pre-main sequence members of associations beyond 500 
pc from the Sun, and potentially out to 2 kpc (for the most 
massive and brightest PMS stars). A major improvement in 
the study of nearby associations will come from the radial 
velocities provided in Gaia DR3. This will be a uniform 
data set that will facilitate studies of the kinematic and dy- 
namic properties of both the clustered and diffuse popula- 
tions of associations. For example, between 40-50% of the 
= 40,000 pre-main-sequence star candidates published by 


(2018) will have radial velocities in Gaia DR3, 


compared to the current value of 2%. 


6.2. Multi-object spectroscopic surveys 


Spectroscopy is important in the study of associations for 
verifying the youth of candidate young stars (Section 2-2). 
providing accurate stellar parameters (allowing stars to be 
placed in the HR diagram and stellar ages estimated), mea- 
suring precise radial velocities (to facilitate full 3D kine- 
matic studies) and abundances (for chemical tagging). 

The next-generation of multi-object spectroscopic in- 
struments will be ideal for studying associations due to their 
combination of wide fields of view (typically several square 
degrees) and large numbers of fibres. This will allow them 
to observe hundreds or thousands of stars at a time spread 
over large areas, thus making them ideally suited to the 
study of associations. 

The first such facility to come online is SDSS-V 
(Kollmeier et al.|2017), whose northern and southern hemi- 
sphere telescopes have already embarked on am ambitious 
survey of the Milky Way and the Local Volume, includ- 
ing many thousands of both high- and low-mass young 
stars. The telescopes operate in both low- (R = 2000) 
and high-resolution ( = 20, 000) in the optical and near- 
infrared, with —500 fibres distributed over 2.8 (northern) 
and 7 (southern) square degree fields of view. 

WEAVE (WHT Enhanced Area Velocity Explorer, [Dal] 
is the new multi-object spectrograph being 
installed on the William Herschel Telescope, with first light 
expected late 2022. WEAVE has a 3 square degree field of 
view and ~ 1000 fibres that can observe in low- (R = 5000) 
or high-resolution (R = 20,000) modes. WEAVE will un- 
dertake 8 pre-planned surveys during its first 5 years of op- 
eration, including surveys targeting dispersed young stars 
and associations in the Galactic Plane, and open clusters. 

In 2023, these instruments will be joined by ESO’s 


4MOST (4-metre Multi Object Spectroscopic Telescope, de] 
spectrograph, which has a 4 square de- 
gree field of view and ~2400 fibres that observe in both 
low- (R = 5000-7500) and high-resolution (R = 20, 000) 
modes synchronously. 4MOST will also operate surveys on 
a 5 year timescale, with multiple surveys of young stars, 
clusters and associations in the first slew of surveys. 


6.3. X-ray telescopes 


The eROSITA X-ray telescope was launched in 2019 and 
is performing a survey of the whole sky over 4 years, in the 
form of 8 successive all-sky maps taken over 6 months. The 
accumulated sensitivity to 0.2-8 keV X-rays will be about 


25 times that of previous all-sky X-ray surveys (Predehl 
let al.2021). This has the potential to identify many can- 


didate young stars across the broad areas covered by asso- 
ciations, by virtue of their magnetic and accretion-related 
activity (see Section[2. I}. 

The first results from eROSITA, based on the first scan 
of the sky, have now been released and a first look at the 
Sco Cen OB association was presented by 
(2021). Already the sensitivity is good enough to detect M- 
type association members at ~ 120 pc with saturated lev- 
els of X-ray emission. Extrapolating to end-of-survey sen- 
sitivities, one can estimate that most low-mass PMS stars 
(0.2 € M/Mọ < 2) will be identified out to ~ 500 pc and 
solar-type PMS stars should be detected well beyond 1 kpc. 
Since X-ray selection is unbiased by kinematic considera- 
tions, it is likely that these data will play a very important 
role in the study of the demography and structure of nearby 
associations over the coming decades. 


7. Summary 


OB associations are important objects to study as they 
represent the transition phase between clustered (or unclus- 
tered) star formation and the population of mature stars or- 
biting in the Milky Way. As such they are important not 
just for understanding the evolution of stellar clustering, but 
also provide samples for studies of stellar evolution (at both 
low and high masses) and the evolution of multiple systems, 
protoplanetary disks, and young planetary systems. 

Identifying members of associations in a reliable and un- 
biased fashion is difficult due to their low densities and 
the large areas of sky they cover. Candidate young stars 
are often identified using a combination of optical and in- 
frared photometry, astrometry, X-ray observations and/or 
variability. Validation of youth can be achieved using spec- 
troscopy, by confirming the effective temperatures of high- 
mass stars or by measuring either the abundance of at- 
mospheric lithium or surface-gravity sensitive spectral fea- 
tures in low-mass stars. Once a reliable sample of young 
stars has been obtained, they can be divided into groups or 
sub-groups using multi-dimensional (spatial and kinematic) 
data and a variety of clustering algorithms. 

Recent work has uncovered significant levels of spatial, 
kinematic and temporal substructure in associations. This 
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substructure has helped resolve the large-scale structure of 
associations, explain kinematic anisotropy and complexity, 
and resolve age spreads commonly observed in associa- 
tions. Evidence for the expansion of associations has been 
mixed, with early studies finding no or limited evidence for 
their expansion, while more recent studies (that break down 
the structure of associations using spatial and kinematic in- 
formation) finding that expansion is a common feature of 
associations, though this expansion can be quite anistropic. 
This in particular has implications for whether kinematic 
ages can be used as a model-free age indicator. The large- 
scale star formation history within (and between) associa- 
tions does not show any obvious patterns that would sug- 
gest a single mechanism responsible for the propagation of 
star formation (e.g., triggering of some sort) and it may be 
that multiple mechanisms are at work at the same time. As 
associations disperse they seem to form extended, filamen- 
tary structures that can survive as coherent moving groups 
for up to ~300 Myr. The origin of their filamentary mor- 
phologies may be a combination of primordial structure and 
the observed anisotropic expansion patterns. 

These observations suggest a picture in which associa- 
tions form as extended, highly substructured and with a low 
average density (albeit with some over-densities). Some of 
the over-densities may be dense enough to form clusters, 
and some of these may survive as long-lived open clus- 
ters, while other expand and disperse, often anisotropically, 
thereby preserving an extended and filamentary structure 
for several hundreds of Myr. 

The low densities at which many stars in associations 
form has implications for the star formation process and 
the formation and evolution of binary and multiple systems, 
protoplanetary disks, and young planetary systems. Full 
observational tests of the predictions made by such implica- 
tions have yet to be carried out, but preliminary studies sug- 
gest some differences in the products of star formation be- 
tween stars in dense clusters and low-density associations. 
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